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Abstract. The propagation of electrons in static and uniform electromagnetic fields 
is a standard topic of classical electrodynamics. The Hamilton function is given by 
a quadratic polynomial in the positions and momenta. The corresponding quantum- 
00 ■ mechanical problem has been analyzed in great detail and the eigenfunctions and 

time evolution operators are well-known. Surprisingly, the energy-dependent coun- 
terpart of the time-evolution operator, the Green function, is not easily accessible. 
However in many situations one is interested in the evolution of a system that started 
with emitted particles that carry a specific energy. In the following we present a 
suitable approach to study this type of matter waves arising from a localized region 
in space. Two applications are discussed, the photodetachment current in external 
] fields and the quantum Hall effect in a fermionic electron gas. 

O ■ 
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o ' 

g ' 1. Introduction 

In quantum mechanics, static and uniform electric and magnetic fields 
are represented by a quadratic Hamiltonian (i.e. a second order poly- 
q , nomial in the canonical coordinates and momenta pi). Moshinsky 

and Winternitz carried out a detailed group-theoretical analysis of 
this class of Hamiltonians and their eigenfunctions [40]. Quadratic 
rN \ Hamiltonians are connected to linear canonical transformations [39] 

and therefore a general phase-space approach gives valuable informa- 
tion for their classification in different dimensions. Nieto used the Moyal 
phase-space representation to develop a general method for deriving the 
time-evolution operator for quadratic Hamiltonians [41]. However, its 
energy-dependent counterpart, the Green function, withstands such a 
systematic analysis and is not available in analytic form for many phys- 
ical relevant potentials. Also other methods, like the Feynman path- 
integral approach, are not capable to derive the exact energy-dependent 
Green function. 

In experiments, the energy of particles is often easier controlled than 
the time of travel. Under these circumstances, the energy-dependent 
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Green function is relevant for the description of the system. Monochro- 
matic particle sources arise in numerous applications of quantum me- 
chanics. In accelerator physics sources located far away from the scat- 
tering region lead to boundary conditions in the form of incoming 
plane waves. In this contribution, we study the behaviour of spatially 
localized electron sources in perpendicular, homogeneous electric and 
magnetic fields. Our discussion is based on the framework of quantum 
source theory, a variant of the scattering formalism that is especially 
suited to describe scattering events restricted to a region of finite vol- 
ume [9, 11, 36, 34]. The idea of quantum sources was first promoted by 
Schwinger [47] but has not found widespread attention. Therefore, we 
briefly introduce the concept and some basic results derived from it, and 
stress its connection to the propagator approach to quantum mechanics 
[18]. In fact, stationary elastic scattering at pointlike sources is fully 
described in terms of the energy Green function. Several fundamental 
properties of the quantum system, like the scattering wave function, 
current density distribution, cross section, and local density of states, 
immediately follow from this functional. Here, we explore in detail 
isotropic point sources in crossed external static fields both in two- and 
three-dimensional configuration space. The results are in agreement 
with experimental findings in a recent photodetachment experiment, 
and offer an alternative viewpoint towards the anomalous Hall effect 
observed in low-dimensional semiconductor devices. 



In preparation for our later discussion, we illustrate the quantum source 
formalism using potential scattering in external static fields as an ex- 
ample. We assume that the potential V(r) represents a localized dis- 
turbance, and that the charged quantum particles otherwise move in 
the external electromagnetic potentials A(r), <I>(r): 



The scattering solutions which are eigenfunctions of the Hamil- 

tonian H with energy E, then usually are decomposed into two parts, 
tp(r) = tpm(r) + tp sc (r), where the incident wave is a solution for the 
unperturbed system Ho: Hoipi n (r) = Eipi n (r), whereas the remainder 
■0sc( r ) represents the scattering wave. By comparison with (1), we find 
that ipsci?) obeys: 



2. Elastic scattering and quantum sources 





[E-H - V(r)} V S c(r) = F(r)^ in (r) . 



(2) 
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Hence, ip S c( r ) is a solution to the inhomogeneous Schrodinger equation 
of the full Hamiltonian H = Hq + V(t), where we denote the right-hand 
side in (2) as the source term <r(r): 

a(r) := V(r) ^(r) . (3) 

Equation (2) suggests the following physical interpretation: The in- 
coming wave ipi n (r), via the perturbation V(r), feeds particles into the 
scattering wave tp sc (v) that is governed by the Hamiltonian H. While 
not commonly seen in standard quantum theory, inhomogeneous par- 
tial differential equations are familiar from other branches of physics, 
the heat conduction equation and Maxwell's equations being examples 
for the introduction of sources. For these problems, a sophisticated 
mathematical framework in the form of Green functions has been devel- 
oped. Accordingly, we introduce the energy-dependent Green function 
G(r,r';E) for the Hamiltonian H defined via [16] 

[E - H - V(r)] G(t, r'; E) = <5(r - r'). (4) 

Formally, the solution to equation (2) is given by a convolution integral 
comprising the source term and the Green function 

^ sc (r) = J d 3 r'G(r,r';E)a(r'). (5) 

We infer that the scattering wave generated by the source <r(r) allows 
for an interpretation as the linear superposition of "fundamental" waves 
G(r,r';E) emitted from point sources C5{v — r') located at r'. 

2.1. Connection to the propagator 

In the continuous spectrum of H, the Green function is not uniquely 
defined. Depending on our choice for G(r,r';E), we obtain a set of 
wave functions tp sc {r) that differ only by eigenfunctions ^hom( r ) of H. 
This ambiguity is resolved by the demand that G(r,r';E) presents a 
retarded solution that enforces outgoing-wave behaviour of the scat- 
tering wave V'sc( r ) at large distances from the source. The representa- 
tion of G(r,r';E) as a Laplace transform of the quantum propagator 
K (r, t|r', to) [18] guarantees the proper choice of boundary conditions 
for the Green function [16]: 

G(r,r';E) = -i lim / dT e ^T/H- v T/H R , Q y (6) 

n 11^0+ Jo 

where K(r, t\r' , to) denotes the coordinate space representation of the 
time evolution operator U(£, to) 

tf(r,t|r',t ) := (r|U(t, t )|r'). (7) 
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Since for a conservative system, U(T, 0) 
formally integrate (6) to obtain: 



exp(— iHT/K) holds, we may 



G(r,r';£) 



lim 



1 



E - H + it] 



(8) 



Therefore, the Green function represents the resolvent assigned to the 
Hamiltonian H in configuration space. According to (8), G(r, r';E) 
indeed acts as an "inverse" to the operator E — H. 

At least in principle, knowledge of the full Green function permits 
the exact evaluation of the scattering wave 4>sc( r ) (5). In general, how- 
ever, G(r, r'; E) is not available in analytic form. In the favourable situ- 
ation that we can find an expression for the Green function Go(r, r'; E) 
associated with the unperturbed Hamiltonian Ho, G(r,r';E) formally 
may be expanded into a series via the Dyson equation: 



1 



1 



E — H E — Hr 



1 + V 



1 



E-U 



(9) 



Replacing G(r,r';E) by Go(r,r'; E), i. e., neglect of the rescattering 
terms that involve the perturbation V, is equivalent to the leading 
order of perturbation theory in the conventional scattering formalism, 
which we will endorse in the following. 

While the quantum propagators K(r, t\r', io) are tabulated for a 
fairly extensive set of potentials [21, 32], few energy Green functions 
are available in closed form for problems in more than one spatial 
dimension. This list includes free particles in two and three dimensions, 
as well as the Coulomb problem [26, 27], uniformly accelerated quantum 
motion [10, 13, 19, 49], the isotropic harmonic oscillator [3], motion in a 
homogeneous magnetic field [15, 20], and parallel electric and magnetic 
fields [17, 35], all in three-dimensional configuration space. 



2.2. Currents generated by quantum sources 

A first obvious quantity of interest are the currents associated with the 
scattering wave ip sc (r) generated by the source <r(r) (5). The current 
density distribution j(r) is defined in the usual fashion via 

j( r) = A^ sc (rrVVsc(r)] - ^^|Vsc(r)| 2 , (10) 
m m 

where A(r) denotes the vector potential, and displays the spatial distri- 
bution of the quanta in the scattering wave, i. e., is directly related to 
the differential cross section of the scattering process. Integration of j(r) 
over a surface enclosing a(r) will yield the total current J(E) emitted by 
the source which, in turn, is a measure of the total scattering rate. For 
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a concise expression, we first note that the inhomogeneous Schrodinger 
equation (2) gives rise to a modified equation of continuity. Instead of 
V • j(r) = 0, valid for a stationary system in the absence of sources, we 
now find: 

V-j(r)=-^[cT(r)*Vsc(r)]. (11) 

Thus, the inhomogeneity cr(r) acts also as a source for the particle 
current j(r). Since the current is conserved outside the source region, 
the surface integral may be replaced by a spatial integration over V • j (r) 
covering the source volume, and upon insertion of (5) for the scattering 
wave in (11), we obtain a bilinear expression for the total flux J(E): 



j(E) = -fa 



J d 3 r J dVcT(r)*G(r,r';£)cT(r') 



(12) 



For pointlike sources a(r) = C5(r — r'), where C is a measure for 
the source strength, the calculation of the scattering currents simpli- 
fies considerably. In this case, the scattering wave is a multiple of the 
Green function, ^ sc (r) = CG(r,r'; E), and the pattern of the current 
distribution follows from (10). Point sources yield a particularly simple 
expression for the total cross section: 

J(r'; E)=- 2 ^f- lim 9 {G(r, r'; E)} . (13) 

h r^r' 

In passing, we remark that for r — ► r', the Green function G(r,r';E) 
diverges in more than one spatial dimension, while its imaginary part 
remains well-defined in the limit and is proportional to the scattering 
rate. The statement (13) is closely related to the optical theorem of 
conventional scattering theory [45]. 

For reference, we list the total currents emitted by a free-particle 
point source of unit strength (C = 1) in one-, two-, and three-dimen- 
sional configuration space. For E > 0, they read: 

J iree \P) ~ ^ . J free W ~ ^ ' J frec W ~ ^3 ' ^ 14 ^ 



where k = \/2mE/h denotes the wave number of the particles. 
2.3. Density of States 

Somewhat surprisingly, the local density of states (LDOS) n(r'; E) 
of the quantum system, i. e., the accumulated density (^(r)! 2 of the 
eigenfunctions of the system with energy E, evaluated at r', is, apart 
from a prefactor, identical to the total current J(r'; E) emitted by a 
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point source located at the same position. Formally, this equivalence is 
established from equation (8) by setting r = r' and using the distribu- 
tion relation [24] $s[(z + i?])" 1 ] = — 7rsgn?7 • 5(z) that holds in the limit 
T] — > 0. Thus, we formally obtain: 

Z[G(r',r';E)] = -ir(r'\d(E-H)\r') . (15) 

The right-hand side of this relation formally contains the spatial repre- 
sentation of the density of states operator S(E — H), and we conclude 
that the LDOS is linked to the imaginary part of the Green function. In 
conjunction with (13), this implies that the density of states is directly 
proportional to the previously defined total current J(r';E): 

J(r';E) = T \C\ 2 n(r';E). (16) 

We note here that the localized eigenstates of H that make up the 
discrete spectrum of the Hamiltonian are irrelevant for the imaginary 
part of the Green function, and thus do not contribute to the current. 
It is the unbounded solutions in the continuous spectrum of H that are 
entirely responsible for the outgoing wave character of G(r,r';E) and 
constitute J(r';E). From (16), we conclude that a non-vanishing den- 
sity of states is therefore directly related to an extended flow pattern in 
position space. We will show examples of this behaviour in Section 5.6. 



2.4. Construction of the Green function 



In the following, we will briefly discuss some techniques that are useful 
in establishing the energy Green function for simple systems. 



Matching of solutions For one-dimensional systems, the inhomoge- 
neous Schrodinger equation (2) reduces to a linear ordinary differential 
equation of second order, and if it exists at all, the Green function 
G(z,z';e) is always available as a product of solutions i/) e ,<(z<) and 
■0 e ,>(.z>) that behave regularly in the sectors z — > ±oo, respectively, 
and are matched at the source position z = z'\ 

G{z, z;e)= - =— . 17 

h 2 W[lpe,<, ll>e,>] 

Here we introduced the symbols z< = min(z,z / ) and z > = max(z,z'), 
and W[tp tt< ,ip t!> ] denotes the Wronskian of the two solutions. The 
basic example for this strategy is the free particle problem in one spatial 
dimension, where for E > (E < 0) ip e ^ < (z < ) and , 4j e:> {z > ) are outgoing 
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(evanescent) waves in either direction: 




exp(-K|z - z'\) (E < 0) 



(18) 



exp(ifc|z- z'\) (E>0). 



Here, k = \/2mE/h, k 



again denote the wave number 



of the particle. In passing, we remark that the few higher-dimensional 
Green functions that can be found in analytical form usually have been 
determined by formal extensions of this technique [3, 10, 26, 49]. 

Eigenfunction expansion In section 2.1, we found a formal position 
space representation for the retarded energy Green function as a spe- 
cial resolvent of the Hamiltonian operator H. Expanding (8) into a 
complete set of eigenstates \ip t ) of H (where H|-0 e ) = e\tp t )), we may 
alternatively express the Green function as a sum over all properly 
normalized eigenfunctions ip e ( r ) — (rlV'e) of the system: 



We will encounter an example of this decomposition in Section 3.3.1. 

Complex convolution We noted before that the quantum propagator 
K(r,t\r',to) is generally more easily available than the Green func- 
tion G(r,r';E) [21, 32]. In part, this situation is the consequence of 
the simple composition properties of K (r , t\r' , to) . Assume that the 
(conservative) Hamiltonian H of the system can be written as the 
sum of commuting, lower-dimensional parts: H = Hi + H2, where 
H4H2 = H2H1. Then, the corresponding evolution operator obeys 



U(T,0) = e iHT / s = e iHlT /V H2T / s = Ui(T, 0)U 2 (T, 0), and thus re- 



duces to a product of its constituents. This property is transferred to 
their spatial representations, the propagators: 



where ri, r 2 are the projections of r onto the subspaces of Hi, H 2 . 

Unfortunately, the simple multiplicative property (20) does not ex- 
tend to the energy domain. We may, however, exploit it to derive a 
corresponding statement for the Green functions G u (r u ,r' I/ ; E). Equa- 
tion (6) shows that the evolution operator U(T, 0) and the resolvent 
operator [E — H]~ , which yields the Green function in configuration 
space (8), are linked through a Laplace transform. Since the image 




(19) 



K(r,t\r',t ) = Ki(n,t|r'i,i )i<: 2 (r 2 ,i|r 2 ,io) 



(20) 



proc_tkramer.tex; 2/02/2008; 4:01; p. 7 



8 



of a product of Laplace transforms is represented by the convolu- 
tion integral of the images of the factors, we obtain [E — H] _1 = 
2^: / dE' [E' — Hi] _1 [i? — E' — H^] -1 , or, in position representation: 

i f°° 

G(r,r';E) = — dE' G 1 (r 1 , r' i; E') G 2 (r 2 , r' 2 ; E - E'). (21) 

27T J-oo 

Due to the generally complicated form of the energy-dependent Green 
function, the practical value of this relation is limited. We will, however, 
present an application in the following section. 

3. Matter waves in crossed electric and magnetic fields 

As our example of interest, we study quantum sources of charged par- 
ticles in an environment of homogeneous, static electric and magnetic 
fields £, 3. The Hamiltonian H in this case may be written as the sum 
of commuting parts H||(r||) and H^(r^) in the sense of Section 2.4 and 
reads: 

(22) 

2 

- qr ± ■ £ ± . (23) 

Here, the subscripts in rii and rj_ denote the directions parallel and 
perpendicular to the magnetic field B, respectively. We chose the elec- 
tromagnetic potentials A(r) = \{S x r) and $(r) = —r£ as particular 
gauge in H = Hii (ri|)+H_i_(rjJ (22), (23), but all observable quantities, 
e. g. the currents j(r) and J{E) (10), (12), are invariant under gauge 
transformations, unlike the propagator and Green function. (We note 
that under a change of gauge field x( r >^)) the source term (3) must be 
correspondingly modified.) 

3.1. The quantum propagator 

According to the composition properties outlined in Section 2.4, the 
propagator for a particle in homogeneous fields £, 3 at arbitrary angle 
may be written as a product (20): 

K £iB (r,t|r',0) = K ll (r ll ,t\r\ ] ,0)K ± (r ± ,t\r' ± ,0) , (24) 

Here, K\\ (ry , t\r'^ , 0) is the propagator for a uniformly accelerated par- 
ticle in one dimension that has been known from the beginnings of 



H l|( r l|) 
Hi(ri) 



— — - r£ 



1 

2m 



II 

-iW 



-(S x r_ 
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quantum mechanics [10, 18, 31, 32]: 



*qi(r||,t|r[|,0) = y 



2iriht 



(25) 



xexp - 



2c2 + 3 



— C — rii ) H cii(ri|+rii) — 

2V 11 H 2 " " I' 24m 



whereas ETj_(rj_, t|r^_, 0), the propagator for a charge moving in two 
dimensions subject to perpendicular electric and magnetic fields, was 
unraveled much later [25, 41, 14]: 



K ± (r L ,t\r' ± ,0) 



ex Pi - 



|B • (r± x r± ) 



(26) 



cot(w L t) [(r ± - r' ± ) 2 - 4v D • (r ± - r' ± )t + v 2 D t 2 



In this expression, we introduced the Larmor frequency ojl and the 
drift velocity V£>: 



wl = qB/(2m) 



v D = (Sx B)/B 2 . 



(27) 



Interestingly, Schwinger derived the relativistic propagator much earlier 
[46], but apparently no simple transition to the non-relativistic case 
exists. 

The equations (25) and (26) reveal a simple symmetry property of 
the propagator under translations of the coordinate origin: 



Kej3{r,t\c',ti) = exp 



B • (r' x r) + r' ■ St 



K £>B (r-r',t\o,0) . 

(28) 

(Alternatively, the symmetry (28) may be viewed as the effect of a 
gauge transform that shifts the origin of the potentials [34].) Of pri- 
mary interest in our study is the corresponding energy Green function 
Gs.b{ t i r 'j E) (4). Its analytical expression is unknown, however. Hence, 
the results displayed in subsequent figures were obtained by numerical 
evaluation of the integral representation (6). Fortunately, the symme- 
try (28) of the propagator, together with (6), permits to predict the 
behaviour of G£ t &(r,r';E) under coordinate transformations: 

G £)B (r,r';E) = exp j^S • (r 7 x r) j G £ , B (r - r', o; £ + gr' • £) . (29) 

This relation immediately extends any result obtained for r' = o to 
general source locations r' ^ o. 
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3.2. Purely magnetic field 

As a simple example, we first inquire into the dynamics of a charge in a 
purely magnetic field {£ = o). In two spatial dimensions, the relevant 
propagator K±(r±, t\o, 0) (26) reduces to: 

^ (r ^ t|O ' 0) = 2vrmsinKt) eXP l^r^ COt(u;Lt) j ' (30) 

As the propagator is periodic in t, we may expand it into a series using 
the generating function of the Laguerre polynomials L/%(z) [1]: 

— exp -H- =]TL n (*)/\ (31) 

L Z 1 J n=0 



UJ L t 



where we set z = e 2luJLt and x = mtOLr^/h. This procedure yields: 

4 2D )(r ± ,t|o,0) = !^L e -^ L rl/2n g Ln / rru^\\ e -i(2n + l) 
71-/1 n=0 \ h J 

(32) 

This form is easily recognized as the decomposition of the time evolu- 
tion operator Uj_(t, 0) into the eigenfunctions of Hj_ (23) populating 
the Landau levels E n = (2n + 1)Huil- 

The two-dimensional energy Green function Gg D ^ (rj_, o; E) follows 
from (32) after the Laplace transform (6) which immediately yields an 
infinite series expression: 

4 2D) (r ± ,o;£) = I^V— V lim V L " ^j/S , 
B v - L ' ' 7 vr/i rj^o^^-^ i (2n + l) + i77 

(33) 

which clearly resembles the formal eigenfunction expansion (19). In 
the important case rj_ — > o, we extract the density of states for a 
two-dimensional gas of charges subject to a magnetic field: 



(34) 

Here, we again made use of the distribution relation Q[(z + ir]) ] = 
— 7rsgn?7 • 5(z) [24] that holds for r/ — > 0. The resulting discrete 5- 
array is indicated in Figure 1. Equation (34) expresses the fact that 
the eigenstates take on only the discrete energy values E n . 

In three spatial dimensions, we have to multiply the propagator 
Kg D \rj_, t\o, 0) with the free-particle propagator in one dimension 
K^\r h t\0,0) = v / m/{2Triht) exp[imr 2 /(2fe)] that follows from (25) 
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Energy [units of huji] Energy [units of Hljl] 

Figure 1. Electronic density of states in two (left panel) and three dimensions (right 
panel) in a purely magnetic field. The dashed line denotes Wigner's threshold law 
(14), valid for free particles. 



once we set £\\ = 0. Its Laplace transform (6) is the free-particle energy 
Green function G^J (ry , 0; E) (18) that we derived in the preceding sec- 

L free 



tion. Similarly, we may transform the product of Kr^\r», t\0, 0) with 



the series expansion of K^ D \r±,t\o, 0) (32) to determine the three- 
dimensional Green function of a charge in a homogeneous magnetic 
field: 

G% D \r,o;E) = ^e- m ^i/2n (3g) 
Tin 

CO 

x^ L n (rnu; L rl/h) G^[r p 0;E- (2n + l)hu> L ] . 

n=0 

(A representation in closed form is stated in [15].) Only the terms 
with positive effective energy (E — (2n + l)ftiVL > 0) contribute to the 
density of states rig (E) = —^{Gg D \o,o;E)], as comparison with 
(18) shows. With the help of (14), (18) we find: 

_gp, m VVge[ E -(2n + l) faL ] 

8 ( ' V2*W „t; JE-(2n + l)h,. L 1 ' 

This superposition of effectively one-dimensional free-particle sources 
is displayed in Figure 1. In passing, we point out that the limes B — > 
in (36) is not well-defined; only after averaging over a small energy 

range 5E, the Wigner free-particle law n^^} (E) = mk/2TT 2 h 2 (14) will 
emerge. Finally, we note that the result (35) extends to the case of 
parallel electric and magnetic fields £ || B, once the free particle Green 
function in the sum is replaced by the one-dimensional Green function 
for a uniformly accelerated particle [17, 35]. 
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3.3. Crossed electric and magnetic fields 

Unlike the case of a purely magnetic or parallel fields, the energy- 
dependent Green function Gs^^r, o; E) for a particle in crossed electric 
and magnetic fields £, 3 is not available in closed form. Thus, in 
this section we limit our considerations to the current emitted by a 
point source, or equivalently, the density of states n(o; E). In a purely 
magnetic field, the degeneracy of the energy spectrum leads to pecu- 
liar shapes of the density of states functionals n&(E) (see Figure 1). 
The presence of an additional perpendicular electric field £ lifts these 
degeneracies and renders a broadened Landau level structure. 

3.3.1. Density of states in two dimensions 

We first examine the two-dimensional case. Here, the Green function 
@e x)g( r -L' °! in perpendicular fields is formally given by the Laplace 
transform (6) of the propagator K±(r±, t\o, 0) (26). Since we are only 
interested in the imaginary part of the Green function at rj_ = o (13), 
(16) we may rewrite this relation and express the density of states 
n f xb(°' as ^ e Fourier transform of the propagator: 



n {2D) 



1 r°° 

(o;E) = —J^dTe iET / h K ± (o,T\o,0) . (37) 

(Formally, the density of states operator 5{E — H) (15) is the Fourier 
transform of the time evolution operator U(T, 0) = exp(— iHT/h), 
and the identity (37) follows in configuration space representation.) 

Alternatively, we may determine n^^(o; E) by direct summation over 
a complete set of eigenstates of Hj_(rj_). We will explore both routes 
below. 

Eigenfunction method A complete set of eigenfunctions for a charge 
in perpendicular fields is conveniently determined in the Landau gauge 
A = (— By, 0, 0) [29]. Here, we assume that the magnetic field points 
into the z direction while the electric field component £± is aligned to 
the y-axis, so the charges drift along the x-axis. The corresponding 
Hamiltonian H',(x, y) is given by: 

H' ± (x, y) = ^- {-ihd x + qByf - ^-8 2 y - q£ ± y . (38) 

The eigenfunctions are products of shifted oscillator functions with a 
plane wave in drift direction, and read properly normalized: 

1/2 



^ (mv D - qBy c 



1 fy-Vc 

; U. 



(39) 
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where 

"-«)=(^b) 1/2 ^ a,2H «<«- < 40 > 

Here, the magnetic length / = \Jh/ (qB) determines the extension of 
the wave function in the direction of €±, while the continuous variable 
y c denotes its centroid. H n (£) is a Hermite polynomial of order n [1]. 
The terms in the corresponding eigenenergy E n (y c ): 

E n (y c ) = (2n + l)hu) L + mv 2 D /2 - q£ ± y c , (41) 

reflect the Landau level, the kinetic energy of the drift motion and the 
potential energy in the electric field, respectively. Summation over all 
eigenstates (39) yields the density of states (37): 



4 2 S(o;£) = f; [°° dy c 5[E - E n (y c )] |Vy, c (o)| 2 (42) 

n=0 J -°° 

n=0 v 

Comparison with (40) shows that the density of states is itself a sum 
over squares of regularly spaced oscillator functions, albeit in energy 
space; their width T and shifts E n are given by: 

T = q£±l = £±^, E n = E-{2n + l)hw L -mv 2 D /2. (43) 

Note that the centers of these oscillator states coincide with the Landau 
levels, apart from a constant shift due to the drift motion. As £± — > 0, 
the width T tends towards zero, and the discrete energy levels familiar 
from a purely magnetic field emerge (34). 

Propagator method To obtain the density of states (42) using the prop- 
agator transform (37), we must first extract the function K±(o,t\o,0) 
from (26): 

K ± (o,t\o,0)= .™ L eX p(^Mcot(^)-l]| . (44) 

This expression formally resembles the propagator in a purely magnetic 
field (30), and the generating function expansion (31) again yields the 
series expansion of K±(o, t\o, 0). As a function of the width parameter 
T defined above (43), it reads: 

^ (o , t|o , 0) ^ e -™| L „(g) 
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77 



(2n + l)Hu> L 



mv 



D 



(45) 



Indeed, the density of states in crossed fields (42) follows after term-by- 
term integration of this sum in (37), as can be shown using the Fourier 
transform [1]: 



die 



-(t-\xf 



[H n (i)] 2 = 2 n n\^L n (2x 2 ) 



(46) 



For a more detailed discussion of this approach, see [37]. While the 
method appears unnecessarily complicated for the determination of 



n 



(2D) 
ex 



g(o; E), the propagator formalism clearly offers an advantage when 

evaluating the Green function Gg^(r±,o; E) for rj_ ^ o, since a single 
numerical integration (6) will suffice to obtain the complete Green 
function. We will show examples below. 



Canonical transformation method The result gained by the two previ- 
ous methods is also consistent with the mapping of the original Hamil- 
tonian in crossed fields in equation (38) to the Hamiltonian of a shifted 
harmonic oscillator. Details of the corresponding canonical transforma- 
tion and its unitary representation are discussed in [33, 34]. 

Properties of the density of states functional The functional form of 
the two-dimensional density of states (42) in crossed fields is a major 
topic in Refs. [34, 37]. Here we content ourselves with a short summary 
of the main features. As noted above, the density of states consists 
of a sum over equally weighted harmonic oscillator eigenstates that 
appear not in configuration space but as functions of the energy E. 
As the eigenfunctions u n {£) (40) form an orthonormal set, the total 
contribution of each sum term, i. e., each Landau level is given by: 

/>>&(*) m 

This result is in accordance with the quantization of the Landau levels 
in a purely magnetic field (34). For each Landau level n, the density of 
states has a Gaussian envelope with width T that is split into n+1 inter- 
vals by the n simple zeroes f nj - (j = 1, . . . , n) of the polynomial H n (£). 
In Fig. 2, we plot the resulting density of states for various electric field 
strengths £±. For small £±, the overlap between adjacent Landau levels 
is negligible, as the DOS drops off exponentially between them. With 
increasing electric field, the Landau levels broaden and finally coalesce. 
We infer from equation (40) that the classical turning point of harmonic 
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Figure 2. Two-dimensional local density of states (LDOS) n(o; E) (in units of 
m/(irh 2 )) and integrated LDOS N(p;E) (in units of qB/(2ivh)) at four different 
electric fields £ = 2000, 4000, 8000, 12000 V/m and for a magnetic field B = 5 T as 
a function of the scaled energy E/(!ujl) according to equation (42). Near the nth 
Landau level at E = (2n + l)!kUL, the DOS renders the probability distribution of 
a one- dimensional harmonic oscillator in the nth eigenstate. 



motion, £* p = y/2n + 1, provides a practical measure for the width of 

the partial density of states n-nExB^)- The populated region in energy 
between adjacent Landau levels n — 1, n is then approximately given 
by the ratio: 

combined half widths of levels r(£* p i + £* p ) , — m£\ , . 

— — ~2V2n - ± . 48 

level spacing 2nwL \j qhB 3 

The overall extension of the modulated Landau levels increases with 
ra 1 / 2 for fixed fields. Note that all features of n nEx&(E), including the 
nodes, scale linearly in width with the electric field £±. While for small 
ratios in (48) individual levels remain well separated, with increasing 
overlap the density of states becomes a smooth function of the energy 
E. This transition is clearly visible in Figure 2. 

3.3.2. Extension to three dimensions 

In three dimensions, besides the Green function G£ t &(r, o; E) even the 
density of states functional ^ (o; E) defies evaluation in closed form. 
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E [units of Hlol] E [units of Hlol] 

Figure 3. Electronic density of states in two (left panel) and three dimensions (right 
panel) in perpendicular electric and magnetic fields (B = 0.5 T, £± = 200 V/m). 



However, simple integral representations are available. Starting from 
the identity (37), we may employ the composition property (24) in 
order to obtain an integral representation: 

4 3 B (P;E) = -L |^dTe ii?T /^ ± (o,r|o,0)K||(0,r|0,0) . (49) 

Alternatively, we may formally perform the integration in (49), which 
leads to a convolution integral of the individual transforms, similar to 
(21). This approach yields a simple composition theorem for the density 
of states: 

/OO 
dE'n^l(o;E')4 D \0;E- E) . (50) 
-OO 1 

While the former expression is better suited for numerical calculations, 
(50) yields more physical insight: Inserting the series expansion (42) 
into (50), we infer that the three-dimensional density of states can again 
be interpreted as a sum over individual Landau levels n, where their ac- 
tual contribution n[ 3 gg(o; E) follows from convolution of the oscillator 
function [u n (E n /T)] 2 (40) with the one-dimensional density of states 
rig (0;E — E). (In the case of perpendicular fields (£y = 0), these 
integrals can be expanded into series of parabolic cylinder functions 
[7], but we will not elaborate this point further.) In Figure 3 we com- 
pare the analytic two-dimensional solution (42) and the corresponding 
three-dimensional density of states (49) in orthogonal fields. The close 
relation between both functionals is clearly displayed, as well as the 
separation of \o;E) into individual Landau levels. 
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3.4. Spin 

A slight complication occurs if the motion of charges with spin, like 
electrons, is considered, since the spin interacts with the magnetic 
field 3. However, for uniform magnetic field, this interaction merely 
causes a constant effective energy shift AE = ±\ghw L if we select 
the magnetic field direction as axis of quantization. Thus, the Green 
functions for each spin component follow from its scalar counterpart 
by adjusting their energy, G+ i (r, r'; E) = G(r,r';E ± AE 1 ). Similarly, 
the spin dependent densities of states become 

n u (E)=n(E±±gfa L S ) , (51) 

and the total density of states including spin can be mapped back to the 
scalar quantity: n^i(E) = n-]{E) -j-n^E). Hence we defer the inclusion 
of spin for the moment. 



4. Application: Photodetachment 

In a photodetachment experiment, electrons are detached from nega- 
tively charged ions A~ due to the interaction with a laser field: 

A" + hv -> A + e~ . (52) 

The detached electron is emitted with a definite energy E given by 
the difference between its binding energy (or affinity) and the pho- 
ton energy. This process allows a description in terms of quantum 
sources. In near-threshold detachment [E — > 0), it is reasonable to 
model the ion as a point source because its size is small compared to the 
de Broglie wavelength of the emitted electron. (Here, we consider only 
the generation of s-waves. For the general case of multipole emission, 
see Ref. [11]-) The photodetachment current in external fields is then 
linked to the relevant energy-dependent Green function: For a laser 
beam illuminating the ions for the duration T, their survival probability 
is given by 

R{E) oc exp[—J(E)T], (53) 

where J{E) denotes the total current defined in equation (13). In 
practice, an external electric field can provide a virtual double-slit 
environment that allows to probe the energy of the emitted electron 
(and thereby the electron affinity of the ion) with extreme accuracy 
[5, 6, 11]. 

The combination of electric and magnetic fields imprints a non- 
trivial structure on the detachment rate and allows to identify features 
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Figure 4- Thermally averaged curves for the total photocurrent as a function of 
the electron energy. Magnetic field: B = 1.07 T, ion mass: m = 32 u. Solid line: 
T = 400 K, dashed line T = 950 K (see also [7], Figure 2). The substructure and 
broadening of the Landau levels due to the perpendicular electric field is visible. 
However, the features are washed out (compared to Figure 3) due to the averaging 
over a wide range of electric field values. 



of the underlying energy-dependent Green function. Unfortunately a 
direct experimental observation of these features is obscured by sev- 
eral effects. Typically, the negative ions are confined in an ion trap, 
where they still have a large kinetic energy. In a thermal ion cloud the 
momentum distribution P(p) is given by Maxwell's expression: 

F(p> =^^P exp (- ^,2/(2mteT, )• (54) 

In an external magnetic field, the charges will experience an electric 
field in their rest frame due to the transformation of the fields [28] that 
accounts for the Lorentz force, 

£l(p) = — pxB. (55) 
m 

This electric field is exactly perpendicular to the momentum and the 
external magnetic field. A stationary source will only emerge if we 
consider the photon-electron interaction in the rest frame of the ion. 
Hence, we employ the three-dimensional Green function for crossed 
electromagnetic fields to describe the photodetachment of moving ions 
in a purely magnetic field. 

The averaging effect of varying electric fields due to the thermal 
motion is displayed in Figure 4. A comparison with the plot for a single 
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Figure 5. Ratio of surviving ions R(AQ) in photodetachment of S~ in an external 
magnetic field as a function of the laser detuning AQ. The solid line is the theoretical 
prediction, the circles represent experimental data [56]. Parameters used: 6 = 1 T, 
ion mass m = 32 amu, T = 2100 K. 

value of the electric field (right panel in Figure 3) shows that the sub- 
structure of the Landau levels changes. Some features are still visible, 
like the division of the first level due to the zero of the first Hermite 
polynomial. Another complication stems from the Zeeman splitting of 
the ionic energy levels in an external magnetic field. Experimentally, 
usually a superposition of many allowed transitions is observed. A re- 
cent experiment is compared to the theory sketched here in Ref. [56] . As 
shown in Figure 5 the agreement is excellent and underlines the validity 
of the quantum source approach. An alternative theoretical description, 
together with earlier experimental data is put forward in [7, 8]. 

5. Application: Quantum Hall effect 

Another application of the Green function in crossed fields is the quan- 
tum Hall effect in a two-dimensional electron gas (2DEG). In Figure 6 
we show the basic geometry of the sample. In the system, a constant 
current is sent along the x-axis of the sample. Perpendicular to the sur- 
face of the electron gas a strong magnetic field is applied. In a classical 
picture, initially electrons entering the sample are deflected to one edge, 
and a potential across the sample builds up until the Lorentz force is 
compensated by the induced electric field. The electrons then drift in 
the crossed fields with the constant velocity vq (27) perpendicular to 
both fields. (Note that this mechanism leads to a loss-free stationary 
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magnetic field 




Figure 6. Schematic picture of a Hall bar. A constant current I x is flowing along 
the x-axis. Perpendicular to the current and an external magnetic field, the Hall 
field is established along the y-axis to counterbalance the deflection of the electrons. 
Experiments record the Hall potential U y and the longitudinal potential U x . 



current in the presence of an electric field, unlike conventional transport 
theory, where the current is limited by inelastic scattering instead.) 
The linear relation between electric field and current density in two 
dimensions is expressed by the conductivity tensor <x: 

( jx \ _ ( &xx °~xy \ ( £x 
V Jy ) V a yx °~yy J \ &y 

Its inverse, the resistivity tensor p, is related to the conductivity via 

°~yx °~xx (K7\ 

Pxy — j T~ i Pxx — j 2~ ■ \?>) 

w xx 1 w xy w xx 1 " xy 

We remark that for a xy ^ 0, vanishing resistivity p xx implies vanishing 
conductivity a xx . In the setup shown in Figure 6, for stationary current 
density j x and transverse electric field £ y , the current j y has to be 
zero. We will now explore different models for the conductivity in the 
two-dimensional Hall effect. 

5.1. Drift transport of electrons 
5.1.1. Classical transport 

In this section we will review the classical transport of electrons in a 
Hall sample. We will take a modest level of scattering of the conduction 
electrons into account. In a simple Drude-like model the dynamics of 
the electrons is governed by the Lorentz force, amended for a term that 
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Figure 7. Classical trajectories from a point source © located at x = y = 0. The 
magnetic field B = 5 T is oriented perpendicular to the plane, the electric field 
£ = 4000 V/m along e y . Left panel: Energy E = hui^, right panel: E = 3huj£. While 
the radius of cyclotron motion varies, the average drift velocity vd — £/B is the 
same in both panels. 



incorporates friction via a relaxation time r: 

dv m , . 

m— =et + evxe v. (58) 

di r 

Under stationary conditions v is constant and together with the cur- 
rent density j = iVev the components of the conductivity tensor in 
equation (56) become 



e 2 Nr 1 / 1 



m 1 + uj^t 2 \ ujct 1 



-UJ C T 



(59) 



where ojq = = eB/m. Inverting this matrix we extract the resis- 
tivity components 

13 Tn 

P*y = N~ e > = j^r T ( 60 ) 

In order to connect this picture to the classical drift of electrons in 
crossed fields we use the relation 

B 

£y = Pxyjx = Jj-3x, jy = 0. (61) 

Solving for v% = j x /(eN) yields v% = vd = £y/B: The current density 
along the x-direction is given by the electron density N, multiplied by 
the drift velocity vd- Some classical electron trajectories are shown in 
Figure 7. As we will see in the next section, the quantum mechanical 
picture radically diverts from these results. 
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Figure 8. Current density from a point source © located at x — y = in crossed 
fields B = 5 T, £ y = 4000 V/m (cf. Figure 7). The electron energy is E = hw L . Left 
panel: The component jx{x,y) at x = —0.1 /um. Right panel: Spatial current flow. 
The arrows indicate the direction of the current and their length is proportional to 

bf /4 



If we assume a field independent carrier density N, equation (60) 
leads to the classical Hall effect: The Hall resistivity p xy is linearly 
dependent on the magnetic field B, and the constant of proportionality 
renders the number of carriers that participate in the transport. 

5.1.2. Quantum mechanical drift 

What do the quantum current and transport look like in our quan- 
tum source model of electronic matter waves? Equations (10) and (13) 
yield quantum mechanical expressions for the currents originating from 
a point source. In Figures 8 and 9 we plot the spatial current dis- 
tribution generated in a magnetic field of B = 5 T and an electric 
field of Sy = 4000 V/m. In the first plot we chose E = huj^ which 
corresponds to the first maximum in the density of states (42) (see 
Figure 3), whereas in the second plot the energy E = 3Hujl is close to a 
minimum in the total current. Some corresponding classical trajectories 
are shown in Figure 7. The quantum mechanical current distribution 
shows some intriguing features: In the vicinity of the source (located at 
the origin), a complicated flow pattern emerges. At some distance from 
the source, the current follows the classical drift direction, but is split 
into two stripes with anti-parallel current vectors. We will discuss the 
implications of these oppositely flowing currents again in Section 5.6. 

The current distribution for stronger Hall fields is depicted in Fig- 
ure 10. Here the drift velocity is vjj = 4000 m/s. Three classical 
trajectories are included in the plot. The nearly circular orbits of Fig- 
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Figure 9. Current density from a point source © located at x = y = in crossed 
fields B = 5 T, £ y = 4000 V/m (cf. Figure 7). The electron energy is E = 3hw L . 
The arrows indicate the direction of the spatial current flow and their length is 
proportional to \j\ l ' 4 . 



ure 7 are distorted to trochoidal shapes [43], and also the quantum 
mechanical current profile changes considerably. 

Total current vs. current density The equation of continuity is valid 
for any surface enclosing the point source, and therefore the spatial 
current density integrated over such a closed surface must yield the 
total current, which by (16) is proportional to the density of states (42). 
Since the total current is available in analytic form, we may use this 
relation to cross-check our numerical evaluation of the spatial current 
density. As a function of E, the functional (42) repeatedly virtually 
drops to zero. This almost vanishing total current does not imply van- 
ishing current density, however, as Figure 9 demonstrates: Circular flow 
patterns lead to a very small net-flow from the source region. Further 
away the contributions of oppositely directed flows cancel each other 
almost perfectly. 



Recovering the drift velocity In order to define an average velocity 
along the drift direction we proceed as follows: For a slice along the y- 
axis at some fixed distance x from the source we calculate the integrated 
density g(x; E) using 

/+oo 
dy |G(r= (x,y),o;E)\ 2 (62) 
-co 
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Figure 10. Current density from a point source © located at s = j/ = in perpen- 
dicular fields B — 2 T, £ y = 8000 V/m. The electron energy is _E = The arrows 
indicate the direction of the current flow and their length is proportional to 
Three classical trajectories of electrons with the same energy are also shown. 

and then define the ratio 

vf{x-E) = J{o-E)/ Q (x-E) (63) 

as the average velocity. This procedure yields values very close to 
the classical drift velocity vo = 800 m/s in both cases illustrated. 
However, J(o; E) and g(x; E) are drastically different for the two en- 
ergies chosen in Figures 8 and 9. (We note that the local velocity 
field v(r) = j(r)/|G(r = (x, y), o; E)\ 2 greatly varies with r. Only the 
integrated quantity reproduces the drift.) Furthermore, it is important 
to realize that the intensity of the current J(o; E) is exponentially 
suppressed at certain energies as shown in Figure 3. This is in sharp 
contrast to the classical picture, where a constant drift transport occurs 
for all energy values of the injected electrons. 

5.2. Fermionic matter waves 

The density of states is a single-particle quantity. In a solid, many 
electrons take part in the conduction process. For a non-interacting 
system, the available single-particle energy levels are occupied accord- 
ing to Fermi-Dirac statistics. Taking spin into account, two electrons 
may share each state. For a system that exhibits a point spectrum of the 
energy levels (e.g. atoms, or a purely magnetic field in two dimensions 
(34)), the resulting electronic configuration is similar to the shell struc- 
ture of atoms. For a continuous spectrum, the Fermi energy controls the 



proc_tkramer.tex; 2/02/2008; 4:01; p. 24 



Propagation in crossed electric and magnetic fields 



25 



integrated carrier density of the system, which at temperature T — ► 
is given by the integrated density of states N(o;Ep)'- 



In crossed electric and magnetic fields, (64) is available in closed form; 
see Ref. [34], Appendix C. Clearly, N(o;Ep) depends on the external 
magnetic and electric fields. 

Here, we should point out that in multi-electron systems Coulomb- 
type interactions will occur. In the following we will ignore these in- 
teractions. The Pauli principle (which requires antisymmetric wave- 
functions for fermions) ensures that two electrons sharing the same 
spin will not be at the same position. Moreover scattering events that 
redistribute the electrons to different energy states take place only if 
the process involves filled initial and empty final states. At T — > 0, 
these states are available only close to the Fermi level of the system. 

5.3. Fermi energy in open and closed system 

Closed system If we treat the two-dimensional electron g ELS clS cl system 
that is closed and decoupled from reservoirs, the number of electrons 
is a fixed quantity. The Fermi energy is determined from the relation 



N = const = iV(o; E F , £ y , B) = / dE n(o; E,£ y ,B). (65) 



Changing the external fields (and therefore the density of states n(o; E)) 
leads to jumps in the Fermi energy as depicted in Figure 11 (top left 
panel): The upper limit of the integral has to be adjusted in order 
to keep the number of carriers N constant. Another implication of 
constant carrier density is a linear relationship between p xy and B 
(60). Thus different mechanisms have to be invoked in order to ex- 
plain the existence of finite Hall "plateaus" of constant conductivity in 
the closed system. Proposals include the formation of one-dimensional 
conduction channels along the edges of the sample, and disorder. In 
a one-dimensional device Landauer quantization gives rise to discrete 
values of the conductivity. Disorder is supposed to lead to localized 
states populated by electrons which do not participate in the transport 
but nevertheless allow to adjust the Fermi energy smoothly [22]. 

Open system In the following we consider the implications of an open 
Hall system, where electrons can enter and leave the system through 
the contacts. In this picture the Fermi energy is fixed, while the number 




(64) 
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Figure 11. Top left: Fluctuation of the Fermi energy as a function of the magnetic 
field for fixed carrier concentration. Top right: Fluctuation of the current carrier 
concentration as a function of the magnetic field for fixed Fermi energy. Lower 
panels: Corresponding resistivity plots (the dashed line denotes the average number 
of carriers) . 



of particles fluctuates around the average free-particle value observed 
for B = 0. In Figure 11 (upper right panel) we plot the oscillations of 

N(o;E F ,£,B)/N { ™\E F ) (66) 

as a function of the magnetic field B. For comparison, we also show 
the resulting Hall resistivity p xy = B/(Ne) in Figure 11 (lower right 
panel) together with its classical counterpart obtained by using N = 

N^\Ef) (14). The following section is devoted to a detailed discus- 
sion of this curve. Here, we merely note that the difference in carrier 
density vanishes at the intersection points of both resistivity curves. 
Otherwise, excess charges will be present whose electrostatic interac- 
tion will lead to potentials that subsequently alter the Fermi energy 
of the system. In the present discussion we will neglect this feedback 
mechanism. 

Experimental evidence for fluctuations Experiments show two types 
of fluctuations in quantum Hall systems as a function of the external 
magnetic field: 
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Figure 12. Schematic sketch of a possible variation of the Hall potential and the 
Fermi energy throughout the sample according to [2]. 



Density fluctuations are directly observed in [44] and fit well into the 
picture of an open system. These results contradict the basic theoretical 
assumptions for quantum Hall systems in Ref. [22, 55], where N = const 
is used to determine the currents. 

According to another experiment [53, 54] the electrostatic potential 
measured atop the two-dimensional quantum Hall system fluctuates 
as a function of the magnetic field. These changes are interpreted by 
the authors as variations of the (local) chemical potential. However, 
according to the basic ideas of most quantum Hall theories [22] , disor- 
der should buffer these oscillations and lead to a smooth variation of 
the Fermi-energy. The implications of these observations are profound, 
since, in the view of a closed system, jumps of the chemical potential 
would prevent the formation of Hall plateaus [22] (see Figure 11). Dis- 
order must be invoked to allow smooth variation of the Fermi energy. 
From the viewpoint of an open system, these fluctuations are correlated 
with the intersections of the classical and quantum mechanical Hall 
resistivities. The excess charges lead also to a varying electrostatic 
potential above the Hall system. In an open system, disorder is not 
an essential ingredient of plateau formation. A careful analysis of these 
intersection points and their matching to the observed fluctuations is 
crucial to investigate this issue. 

5.4. Fermi energy and Hall potential variations 

Up to now we considered only the injection of electrons from a single 
point source. We will now extend the formalism to cover a continuous 
"wire" of point sources along the current injecting contacts. Also, we 
will introduce the possibility of a local variation of the Fermi energy 
.Ep(y)and the Hall field £ y (y), as sketched in Figure 12. The current of 
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a macroscopic device with width W is given by integrating the current 
density over the width of the device: 



For a given form of the Hall potential and the Fermi energy we are now 
in a position to calculate the resistivity p. 

5.5. Calculation of the Hall resistivity and the current 



For the realistic calculation of a Hall resistivity curve we have to incor- 
porate some material parameters in the theory: 

1. Electrons in solids are characterized by an effective mass m*. 

2. Similarly the magnetic (/-factor of the electron depends on the 
material and possibly on the magnetic field B. 

3. In some materials additional degeneracies appear (e.g. the "valley 
splitting" in silicon). 

4. All observations are made at a finite temperature T. 

5. The electric field and Fermi energy may vary along the direction of 
the Hall field £ y . 

6. In a multi-electron system Coulomb interactions between the elec- 
trons and the positive background charges occur. 

7. In a non-perfect sample disorder and electron-phonon interactions 
are present. 

5.5.1. A new expression for the Hall conductivity 

The combination of the quantum source model with the Pauli principle 
allows us to obtain a purely quantum mechanical expression for the cur- 
rent along the drift direction. Working in the eigenfunction expansion 
for the Green function (see Section 3.3.1), each eigenstate supports the 
current 



£.y c(r) = _Q { ^ yc{rTdx ^ yc{r)} + f-£|^ n)Wc (r)| 2 . (68) 




(67) 



FLOW 



m 



m 



The properly weighted current is given by 




n=0 



(69) 
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A short calculation yields 

CO . 

j x (r;E) = v D J2 / dy c S[E - E n (y c )] \tp n ,y c {r) 



n=0 

eB c 



+— E / d Vc 6 l E ~ E niVc)\ (Vc - y) IV>^ c (r)| 2 . (70) 
m ^o J 

We already evaluated the first term in (42). A (macroscopic) conduc- 
tivity is obtained by integrating over y. Since the second integral runs 
over a function antisymmetric in {y — y c ), the second term vanishes. 
Inserting vp, = E jB gives 



Without integration, the second term is responsible for the complicated 
flow pattern seen in the current pictures analyzed in Section 5.1.2. In a 
many electron system with Fermi energy Ep the conductivity in Ohms 
law j = <r • E becomes 

° xy = I IZ dE n ^ (r; E) = B~ N(r; Ef) ' (72) 
a xx = 0. (73) 

We emphasize that this expression couples the specific form of the 
density of states in crossed electromagnetic fields with the drift velocity. 
It is not possible to separate the quantity iV(r; Ep) from the drift 
velocity and introduce it as an independent classical parameter. 

5.5.2. A simple model for the quantum Hall effect including scattering 
The previous model is not complete, because the longitudinal resis- 
tance is always zero. Experiments show a non-vanishing a xx if Ep 
coincidences with a Landau-level. A natural extension of the model 
is the incorporation of scattering. A simple, yet instructive model for 
the Hall effect that incorporates effects 1-4 is presented in [34, 37]. We 
start from a conductivity tensor similar to the classical expression (59). 
However, we take into account the external field and energy dependence 
of all quantities: 

e 2 n(o;E)r(E) 1 / 1 -ujct(E)\ 

[ >~~ m 1 + c4t(£) 2 \ujot{E) 1 j' [<V 

We will not consider a locally varying electric field. The discussion 
also assumes that the electrons are injected at a point-contact located 
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at o. We can lift this restriction by introducing a position dependent 
Fermi-level [23] in the system 

E F (r)=E F (o) + er-£. (75) 

Usually, unequal Fermi levels result in a current. In the Hall geometry 
this current along the electric field is absent as the electrons only drift 
perpendicular to the magnetic and electric fields. From the symmetry 
relation (29) we obtain a translational invariance in the sense that 

n(r;E F (r))=n(o;E F (o)). (76) 

Therefore we can regard the local density of states as the global density 
of states in the system. A more sophisticated model, which incorpo- 
rates a (slow) variation of the electric fields and the Fermi energy, is 
sketched in Section 5.6. In the Appendix we show, that the Lorentz- 
force model for the conductivity may be replaced by the expression for 
the probability current defined in (13). 

For T — > 0, the total conductivity is obtained by integrating over 
the occupied energy range 

/•Ep 

c = / dEa(E). (77) 



For strong magnetic fields, the energy-dependent relaxation time t(E) 
satisfies [ujc t (E)] 2 3> 1, except in the vicinity of E F , and the transver- 
sal component a xy thus mirrors the integrated density of states iV(o; E F ) 
(64) 

^=irj E i + t'r%)\-i = e B N{ °-' EF) <78) 

[cf. (60)]. The last expression is already known from our first model: 
The quantization of the plateaus in the Hall effect does not depend on 
the scattering. 

The longitudinal component is more difficult to evaluate since it in- 
volves assumptions about the scattering events. If we assume that only 
electrons with energies close to the Fermi energy contribute significantly 
to a xx , we obtain: 



a xx (E F ) = dE n(o; E) « Dn(o; E F ). (79) 

Here D denotes some constant that may depend on the material pa- 
rameters and the fields. For finite temperature, significant scattering 
may also take place in an energy range of several ksT around E F . The 
value of a xx {T) then follows after suitable averaging [52] 

a xx (T) = JdE (~ df f E T) ) (80) 
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Figure 13. Quantum Hall effect for a constant magnetic field B = f9 T. Effective 
mass m* = 0.2m e , T = f.5 K. Notice that any substructure of the Landau lev- 
els is washed out by thermal averaging. The dashed line shows the classical Hall 
line. The lower dash-dotted line is proportional to the longitudinal resistance p xx . 
Corresponding experimental data is shown in [51]. 

where f(E,T) denotes the Fermi-Dirac distribution 



The Hall resistivity for a fixed magnetic field In early experiments 
on the quantum Hall effect, the resistivity was measured under the 
condition of a fixed current along the x-axis (I x = const) and a fixed 
magnetic field [51]. By varying the gate voltage in the experimental Si- 
MOSFET system, the Fermi energy is adjusted. We will assume a linear 
relationship between the gate voltage and the Fermi energy. Then, the 
value of the Hall field £ y is a solution of the implicit equation 

£y = Pxy (E F ,£y,B) j x (82) 

for given Ep,B,j x . Here, the resistivity p xy is related to the conduc- 
tivity components (78), (80) via (57). For the interpretation of data, 
we have to include an additional degeneracy besides spin that occurs 
in silicon, the "valley splitting" which effectively doubles each level, 
leading to a total of four repetitions of each Landau level. Introducing 
the additional valley quantum number v = ±^ and the corresponding 
energy shift E v [50], the density of states given in equation (51) becomes 

«T>l,valiey(o; E) = n (o; E ± gftoj L /2 ± vE v ) . (83) 
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Magnetic field B [Tesla] 

Figure 14- Quantum Hall effect at strong magnetic fields (B > 1 Tesla) for a 
non-interacting two-dimensional electron gas. The plot shows the Hall resistance p xy 
and longitudinal resistance p xx as a function of the magnetic field B for fixed Fermi 
energy (Ef = 0.868 meV). Effective mass m* = m e , effective g- factor g — i, current 
density j x = 0.2 A/m, t(E f ) = 10~ n s, T = 0.1 K. The dashed line represents the 
classical Hall resistance p xy with a constant level density. Experimental results are 
shown in [42]. Note that j x is chosen fairly large in order to show the substructure 
of the Landau levels. 

Temperature dependence If the ratio fcgT '/ '(huii) becomes close to 
unity, the Hall plateaus disappear, since o xx (T) as given by equa- 
tion (80) is no longer approaching zero between two Landau levels. 
However, the introduction of an effective mass m* can lead to large 
modifications of uj* l = eB /(2m*). Since the width ksT is independent 
of material parameters, it can be used as an independent energy scale 
to access the values of m* ,g* . In Figure 13 we assume an effective mass 
m* /m e = 0.2. A higher effective mass would be inconsistent with the 
reported temperature of T = 1.5 K, since a smaller energy range fvjo* L 
cannot accommodate four separated peaks of individual width ksT. 
Thus, the temperature dependence of p xy may be used to determine 
some of the material parameters of the system. 

The Hall resistivity as a function of the magnetic field Nowadays, 
GaAs-heterostructures are commonly used to provide the two-dimen- 
sional electron gas for the quantum Hall effect. Advantages are cleaner 
samples with very high mobilities and the absence of valley splittings. 
However, in these samples the number of electrons for B = is vir- 
tually constant and largely independent of the gate voltage. Therefore 
the magnetic field is varied while keeping the gate voltage fixed. In 
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Figure 14 we show a typical plot of the resulting resistivity. As men- 
tioned before, the intersection points of the averaged or classical Hall 
resistivity with the quantum mechanical result for p xy deliver impor- 
tant information: They allow to analyze the spin splitting and other 
parameters of the system. 

5.5.3. Fractional effects 

As shown before, the electric field leads to additional zeroes in the 
density of states and consequently to subdivided Landau levels. Their 
fractional values of the filling factor are analyzed in [37] and displayed 
in Figure 14. Experimentally, the appearance of the fractional filling 
factors arising from the electric field might be difficult to detect, since 
uniform and high current densities are required at very low tempera- 
tures. Furthermore, the fractional quantum Hall effect may overshadow 
the single-particle structure if Coulomb interactions dominate the Hall 
field contribution. Experiments show plateaus in p xy for simple frac- 
tions of the filling factor in the first (n = 0) Landau level. In the pre- 
sented model such features are not explained. The effect is attributed to 
collective modes of the system that are caused by interactions between 
the electrons [55]. We explicitly did not include Coulomb interactions 
in our model system. 

The best candidates for the detection of field-induced fractional fill- 
ing factors are exactly half-filled odd Landau levels, since they are left 
largely unaffected by the averaging caused by a non-uniform electric 
field. (A similar effect prevails in three dimensions, where the dip 
associated with the second Landau level is still visible in Figure 4 
despite extensive averaging over different field strengths.) In contrast, 
according to the standard theory of the fractional quantum Hall effect 
half-filled levels do not induce plateau formation. 

5.5.4. Hall-field dependence of the plateau width 

Another signature of the presence of an electric field dependent broad- 
ening of Landau levels is the breakdown behavior for high currents. 
Kawaji and co-workers conducted extensive experiments on the charac- 
teristics of the breakdown and find experimentally a dependency that is 
exactly the same as the one obtained in equation (48): The width of the 
plateaus decreases linearly with increasing current. The experimental 
observation of this behaviour is reported in [4, 30, 48]. 

5.6. Current distribution 

The form of the Hall potential is actually experimentally accessible [2] , 
and a schematic result is sketched in Figure 12. We should note that 
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Figure 15. Quantum mechanical current density from an extended "wire" (shaded 
region). The electron emitting region extends from y = to 0.2 /im. Left panel: The 
component j x (x,y) at x = —0.1 /an. Right panel: Spatial current flow. The arrows 
indicate the direction of the current and their length is proportional to |i| 1//4 . In 
this example, an effective edge current forms with unequal magnitudes at each edge. 
The bulk also carries a constant current flow. 



the theory presented here could be applied for any given variation of 
the Hall potential and the Fermi energy. A starting point for models 
that provide this input could be the self-consistent potentials obtained 
in Refs. [12, 38]. For the sake of simplicity, let us discuss here only a 
straightforward extension in which we treat the emission of independent 
electrons along a constant Hall field. We already calculated the current 
distribution for a point source. For the "wire" described above, we 
obtain the global current profile by summing over the current contri- 
butions of the point sources. Here, we will use eq. (75) and assume that 
the Hall field is constant across the probe. 
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Figure 15 displays the resulting flow pattern. The previously de- 
scribed oppositely flowing currents are shifted to the edges, while in the 
bulk a uniform current emerges. In this way effective "edge" currents 
are established in a model of an electron emitting contact of finite 
width. We should note that the magnitudes of the oppositely directed 
edge currents differ, as already seen for point sources. 

6. Conclusions 

The quantum source formalism provides an excellent basis for the anal- 
ysis of the propagation of matter waves in external fields. While some 
classical properties of the motion of particles prevail in the quantum 
mechanical C£LS6 ? EL smooth transition from quantum to classical me- 
chanics is generally not observed. The Landau quantization due to the 
magnetic field and the combination with an electric field have profound 
implications for the spatial current distribution and intensity. 

Present-day nanotechnological devices, operated at very low tem- 
peratures, can actually provide experimental data for the propagation 
of electronic matter waves and allow to test theoretical predictions. We 
analyzed a simple model of the quantum Hall effect as one example. 
Even this non-interacting electron picture already shows a wealth of 
interesting features and gives access to important parameters of the 
system. 
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